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Abstract 

Tidal streams of the Sagittarius dwarf spheroidal galaxy (Sgr) may be showering dark 
matter onto the solar system and contributing ^(0.3-23)% of the local density of our 
Galactic Halo. If the Sagittarius galaxy contains WIMP dark matter, the extra contribution 
from the stream gives rise to a step-like feature in the energy recoil spectrum in direct dark 
matter detection. For our best estimate of stream velocity (300 km/sec) and direction (the 
plane containing the Sgr dwarf and its debris), the count rate is maximum on June 28 and 
minimum on December 27 (for most recoil energies) , and the location of the step oscillates 
yearly with a phase opposite to that of the count rate. In the CDMS experiment, for 60 
GeV WIMPs, the location of the step oscillates between 35 and 42 keV, and for the most 
favorable stream density, the stream should be detectable at the llcr level in four years of 
data with 10 keV energy bins. Planned large detectors like XENON, CryoArray and the 
directional detector DRIFT may also be able to identify the Sgr stream. 



1 Introduction 

The dark halo of our galaxy may consist of WIMPs (weakly interacting massive particles). 
Direct detection experiments attempt to observe the nuclear recoil caused by these dark matter 
particles interacting with the nuclei of the detector. Currently many such detectors are in 
operation worldwide, including CDMS, COSME, CRESST, DAMA, DRIFT, Edelweiss, Elegant 
V, HDMS, IGEX, LiF/Tokyo, LosAlamos, MiBeta, ROSEBUD, Saclay, UKDMC, and Zeplin. 
These detectors may be able to observe not only the general halo component but also streams 
of dark matter moving through our solar system. 
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Recently, the Sloan Digital Sky Survey (SDSS) and the Two Micron All Sky Survey (2MASS) 
surveys have traced the tidal stream of the Sagittarius dwarf galaxy (Sgr) more than 

360° around the sky. It is believed that this tidal stream is currently moving through the solar 
neighborhood, and that it contains dark matter in addition to stars. In this paper we investigate 
the possibility of detecting WIMP dark matter in this stream via direct detection experiments. 

The Sagittarius galaxy, a dwarf spheroidal galaxy of roughly 10 9 M Q , is a satellite of our own 
much larger Milky Way Galaxy, located inside the Milky Way, ~15 kpc behind the Galactic 
Center and ~6 kpc below the Galactic Plane 3. . There are two streams of matter that extend 
outwards from the main body of the Sgr galaxy. These streams, known as the leading and 
trailing tidal tails, are made of matter tidally pulled away from the Sgr galaxy. The leading tail 
may be showering matter down upon the solar system [2]. The flow is in the general direction 
orthogonal to the Galactic plane and has a speed of roughly 300 km/s (see discussion in §2 
below) . This speed is comparable to that of the relative speed of the Sun and the WIMPs in the 
general dark halo. Hence one can hope to detect the stream in direct detection experiments. 

The detectability depends on the density of dark matter in the stream. The mass-to-light 
ratio M/L in the stream is unknown, but is plausibly at least as large as that in the Sgr main 
body; in fact, the M/L in the stream may be significantly larger because the dark matter on 
the outskirts of the main body would be tidally stripped before the (more centrally located) 
stars. Various determinations of the M/L for the Sgr main body give values in the range 25 to 
100 (see the discussion in [2] and references therein). With our assumptions, we find a density 
of dark matter in the stream in the range 0.3 to 23% of the local dark halo density. If the dwarf 
does not contain dark matter or the stream does not pass through the solar neighborhood, then 
the contribution to the local dark matter density is zero. 

This additional flux of WIMPs from the stream shows up as a 0.3-23% increase in the rate of 
nuclear recoils at energies below a characteristic energy E c , the highest energy that WIMPs in 
the stream can impart to a target nucleus. Hence, there is a step in the energy recoil spectrum; 
the count rate in the detector is enhanced at low energies, but then returns to the normal value 
(due to Galactic Halo WIMPs) at all energies above the critical energy E c . This feature can be 
observed as a sharp decrease in the count rate above a characteristic energy that depends on 
the mass of the target nucleus, the mass of the WIMP, and the speed of the stream relative to 
the detector. 

It is well known that the count rate in WIMP detectors will experience an annual modulation 
|S] as a result of the motion of the Earth around the Sun: the relative velocity of the detector 
with respect to the WIMPs depends on the time of year. We find that the count rate of WIMPs 
in the stream is a maximum on June 28 1 and a minimum on December 27. 2 Note that this 
modulation has a different phase than does the ordinary annual modulation of WIMPs in the 
Galactic halo, which peaks on June 2 for high energy recoils (and December 20 for low energy 
thresholds [SI |5J |U| ) ■ Similarly, we find that the value of the characteristic energy E c of the 
step in the energy recoil spectrum oscillates in time with a period of a year. The step location 
experiences an annual modulation exactly 180° out of phase with that of the stream count rate. 
To give a number, for a 60 GeV WIMP detected in a Ge-73 detector, the characteristic energy 
for the Sgr stream varies between 25 and 35 keV. As we will show, this step feature in the energy 
recoil spectrum should be detectable at the three a level in the CDMS experiment using two 

J We suggest that readers prepare their WIMP shields. 

2 For energies near E c , there is a small range of energy recoils in which the count signal has the same phase 
as the characteristic energy and so peaks in December. 
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kilograms of Ge and taking four years of data with 10 keV energy bins. 

In 1994 and again in 2001, the HEAT collaboration |10l detected an excess of cosmic ray 
positrons. Baltz et al. ^1^1 explored the possibility that this excess is due to the annihilation 
of neutralino dark matter in the Galactic halo. They found that neutralino annihilation can 
produce enough positrons to make up the measured excess in the Minimal Supersymmetric 
Standard Model provided there is an additional enhancement to the signal due to dark matter 
clumps in the Galactic halo 3 . Baltz et al. found that, under current constraints on the MSSM, 
the "boost factor" needs to be at least 30 in the signal. The Sagittarius stream provides at 
most a 23% enhancement of the signal, and hence is unlikely to provide an explanation of 
the boost factor required for neutralino annihilation to explain the HEAT data. However, the 
uncertainties in the propagation of cosmic rays in the galaxy are large and it may be possible 
that a smaller boost factor may suffice, so further studies are warranted. 

Detection of this dark matter stream would be important also for our understanding of 
galaxy formation. The traditional assumption that the halo population of stars is old and 
well mixed has been seriously challenged by recent observations that the outer parts of the 
stellar halo contain tidal debris that is merging with the Milky Way even at the current epoch. 
Recent N-body simulations of galaxy formation suggest that there may be dwarf satellites of 
the Milky Way that contain only dark matter and no stars. If dark matter can be detected in 
the Sagittarius stream, then we could use dark matter detectors on Earth to look for remnants 
of these dark companions of our Milky Way. 

2 Local dark matter from the leading tail of Sagittarius 

The factors that affect dark matter direct detection experiments are the local dark matter 
density and velocity vector of material in the leading tail of the Sgr dwarf galaxy that passes 
through the solar neighborhood. Below, we estimate that the (local) Sgr leading tail is moving 
with a velocity of 300±90 km/s in the direction 0.233y — 0.970Z, where Z is the north Galactic 
pole, and Y is in the direction of the local disk stars' motion (I = 90°). Equivalently, the local 
motion of the Sgr leading tail points in the direction (I, b) = (90°, —76°), where / is the Galactic 
longitude and b is the Galactic latitude. All velocities are measured with respect to the rest 
frame of the center of mass of the Milky Way galaxy. We will also find an estimated range for 
the local density of dark matter of: 

PtaU,© = [2, 190] x 10 4 M Q /kpc 3 = [0.001, 0.07] GeV/cm 3 , (1) 

which corresponds to (0.3-23)% of the local density of the isothermal Galactic halo. 

The density calculation is highly uncertain, since (1) the stellar density is not well measured 
anywhere along the stream and has not been directly detected in the solar neighborhood, and 
(2) M/L is highly uncertain in the Sagittarius dwarf itself and is unmeasured in the tidal 
stream. Models of the tidal disruption have not been developed or constrained by observations 
to the point that they usefully predict the local properties of coherent dark matter in the solar 
neighborhood. We will derive our estimates from empirical observations of the old populations 
of stars in the tidal stream. Blue horizontal branch stars and main sequence turnoff stars are 

3 An alternate explanation within the MSSM, namely neutralinos that were produced non-thermally in the 
early Universe, was suggested by Kane et al. 1141 1151 . 
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old enough that their populations are not likely to have significantly changed in the time since 
these stars were stripped from the dwarf galaxy. The important measurable quantities are (1) 
the number of stars observed in some portion of the tidal stream, (2) the volume of that portion 
of the tidal stream, (3) the relationship between number of observed stars of a particular type 
and total luminosity in stars, and (4) M/L in the tidal stream. 

2.1 Velocity in the solar neighborhood 

The Sgr dwarf galaxy is currently about 16 kpc from the Galactic center, near its closest 
approach to the Galactic center (perigalacticon) . As seen from the Sun, Sgr has a radial velocity 
of 140 ± 2 km/s and a tangential velocity of 250 ± 90 km/s |T5]- Correcting for the solar motion 
as given in §3.2, we derive a current velocity of 300 ± 90 km/s for the dwarf, in the reference 
frame of the Galaxy center. 

Stars in the leading tidal tail, which extends to near the solar position, were stripped from 
the side of the dwarf that is closer to the center of the Milky Way. At the time that the stars 
are stripped from the gravitational pull of the dwarf galaxy and released into the potential of 
the Milky Way, they have about the same orbital velocity as the dwarf galaxy itself (we assume 
that the dwarf galaxy stars have no bulk rotation) . Due to their closer proximity to the Galactic 
center, they have a lower gravitational potential energy, and therefore shorter orbital periods - 
thus creating the leading tail. 

Stars in the trailing tail were stripped at larger distances from the center of the Milky Way. 
The higher orbital energy of these stars produces larger distances from the Galactic center and 
longer orbital periods. 

The trailing tidal tail has been traced in Two Micron All Sky Survey (2MASS; [T7j) M giant 
stars through nearly 180° of sky south of the Galactic plane, to a maximum distance of ~50 kpc 
from the Galactic center. The tidal tails lie on a plane that passes through the Galactic center. 
The Sun lies within a kiloparsec of that plane, and near the path of the leading Sgr debris [2] ■ 

The stars in the tidal tails must conserve both energy and angular momentum, in the absence 
of dynamical friction. Since the stars in the trailing tail of Sgr have an apogalacticon of 90 kpc 
or more, while the stars in the leading tail have an apogalacticon of only about 45 kpc p^, the 
trailing tail stars must have significantly higher energies than the stars in the leading tail. If 
most of the stars are stripped near the perigalacticon of the Sagittarius dwarf, at 16 kpc from 
the center of the Milky Way, then the leading and trailing debris must have been stripped from 
parts of the Sagittarius dwarf that differ in Galactocentric distance by at least 10 kpc J5| • It 
is sensible, therefore, that the perigalacticon of the stars in the leading tail could be 11 kpc or 
less. 

Since the mass of the Milky Way is distributed throughout its volume and is not concentrated 
at the center, stars do not travel in closed orbits, but rather in rosette patterns. The angle 
between successive perigalacticon passages is less than 360°. The Sun is near the first expected 
perigalacticon passage of the leading tail, at about the right distance. Since the tidal stream is 
approximately 6 kpc across, it is reasonable to suspect that the tidal steam passes through the 
solar position. Majewski et al. 2 show that stellar debris appears to be raining down on the 
solar position from approximately the direction of the north Galactic pole. 

Without significant dissipative forces, the average speed of the tidal debris in the solar 
neighborhood (which is at perigalacticon) should be the same as the current speed of the Sgr 
dwarf (which is also at perigalacticon). The debris should travel approximately in the plane of 
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the Sgr dwarf tidal debris. If the debris is near perigalacticon, then its motion is approximately 
perpendicular to the line from the Sun to the Galactic center. From these three conditions, 
we derive that the debris velocity is in the direction 0.233F — 0.970Z, where Z is the north 
Galactic pole, and Y is in the direction of the local disk stars' motion (I = 90°). In terms of 
galactic latitude b and galactic longitude /, this direction is equivalent to the direction (l,b) — 
(90°, -76°). 



2.2 Simple estimate of local dark matter density 

In this section we use simple parameters of the Sgr dwarf tidal system to roughly estimate the 
Sgr dwarf contribution to the local dark matter density. In the next section we will make more 
careful estimates of the parameters. 

Published estimates of the original total mass of the Sagittarius dwarf galaxy are in the 
range 10 9 to 10 11 M Q . The lower number is within a factor of two of current estimates of the 
mass of the Sgr dwarf, and the larger number is derived from the need to keep the dwarf intact 
for a Hubble time The width of the tidal stream in 2MASS M stars is estimated to be 4-8 
kpc 0. Below, we will calculate a FWHM of 6 kpc from the data from the SDSS 20 . This 
width is of the order of the size of the dwarf galaxy. 

Estimates of the current mass of the Sgr dwarf galaxy are in the range 

M Sgr ~ [5,20] x 10 8 M Q , (2) 

possibly several orders of magnitude lower than the original mass of the Sgr dwarf at its time of 
formation. The lower mass is the mass within the Sgr dwarf tidal radius calculated by Majewski 
et al. from the dispersion of the radial velocities of Sgr dwarf stars. Comparing their mass to 
the luminosity of the Sgr, they derive a mass-to-light ratio M/L of 25-29 in solar units. Other 
authors have argued for M/L as high as 100 16 (four times as high) based on the need to keep 
the Sgr dwarf intact for a Hubble time. This produces the upper end of the current mass range. 
Given that the M/L ratio is expected to be 25 or higher, calculations of the mass density in 
the Sgr dwarf and tidal tails are essentially equal to the dark matter density. 

A typical disruption model |23 assumes the current mass of the Sgr dwarf is about 10 9 M©, 
and that the dwarf lost about 2-3 times that mass in the last gigayear (somewhat longer than 
an orbit). If this is spread out into a tidal stream 6 kpc in diameter and 250 kpc in length (this 
is the length of the observed tidal stream in stars), then an estimate of the local dark matter 
density is: 

""■■^ (6^g^) = 3Xl0 ' M ^'- 
This is one percent of the estimated local dark matter density. 



2.3 Detailed calculation of the stream density 

We now explore the range of dark matter local dark matter densities that are consistent with 
our knowledge of the Sagittarius system, assuming that the leading tidal tail passes through the 
solar neighborhood. 

In order to calculate the local dark matter density from data, we would like to measure the 
local mass density of stars from the stream, and multiply by M/L calculated for our position on 
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the leading tidal tail. Neither the local density of Sgr stream stars, nor the local mass-to-light 
ratio in the stream, is known. Instead, we can calculate the density of F/G turnoff stars at one 
location along the trailing tidal tail, tabulate the blue horizontal branch star density at several 
positions along the leading and trailing tidal streams, and estimate M/L from the Sgr dwarf 
itself. Using these quantities, we estimate the dark matter density from the Sgr tidal stream in 
the solar neighborhood. 

Specifically, we wish to estimate ptaii,©, the dark matter density of the tail in the neighbor- 
hood of the Sun, as the product of three quantities: 



/ number of stars \ / luminosity \ / mass \ 
Ptaii,© ^ volume J \ number of stars / \ luminosity J 



(4) 



Ideally, one would like to measure all of these quantities in the tidal stream at the solar position. 
Because none of these quantities are known at the solar position, each term will be estimated 
from other measurable portions of the Sgr tidal stream. 

To estimate the first term on the right hand side, we will measure the density of F/G turnoff 
stars in a portion of the trailing tidal tail 50 kpc from the Sgr dwarf (stripe 82 in the SDSS). 
The density in this location is the best measured anywhere in the tidal tails. To obtain an 
estimate of the density at our location, we extrapolate from measurements of Blue Horizontal 
Branch (BHB) stars at three places along the tidal tails. We discover that the density of the 
stars in a piece of the leading tail 50 kpc from the Sgr dwarf (stripe 10 in the SDSS) is nearly 
a factor of two higher, giving us a larger range of possible dark matter densities at the solar 
position: Pf/g.q ~ PF/G,82( p°™'° )■ The second term on the right hand side of Eq.Q is then 

The third term is (M/L) ta ii- Multiplying together these three factors, we find 
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PF/G,82 

Pbhb,© 

PBHB,82 

M Sgr 



-^F/G,Sgr 
M/itail 
Af/Lggr 



the number of Sgr F/G stars per kpc 3 in stripe 82, 

the density of Sgr BHB stars at the solar position, 

the density of Sgr BHB stars in stripe 82, 

the current mass of the Sagittarius dwarf, 

the number of F/G stars in the Sagittarius dwarf, 

the mass-to-light ratio in the tidal tail, 

the mass-to-light ratio in the Sgr dwarf. 



Each of these terms is discussed in further detail in the following sections. 



2.3.1 Estimating /?f/g,82 

First, we measure the stellar density of F/G turnoff stars in stripe 82, (l,b) = (167,-54), 
about 50 kpc from the Sgr dwarf along the trailing tidal tail. Figure 22 of Newberg et al. [2U| 
shows the counts of F/G turnoff stars along the Celestial Equator. From that figure, one can 
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50 kpc, leading 
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0.9848 
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0.67 


150 kpc, trailing 


82 


35 


29 


0.9848 


0.5071 


0.60 


50 kpc, trailing 



Table 1: Calculation of relative A star density at three places along the Sgr dwarf tidal tail. 
Stripes 29 and 82 cross the trailing tidal stream, and stripe 10 crosses the leading tidal tail. 
Stripes 29 and 10 are near perigalacticon for the trailing and leading tidal tails, respectively. 
The last column gives the adjusted relative A star density for each piece. 



estimate that in the Sgr stream overdensity labeled S167-54-21.5 (in red), there are 1000 stars 
per half-degree-wide bin at the peak. It is difficult to estimate the background star counts in 
that region of the sky, but it is between zero and 400 stars per half degree bin. The total width 
of the (approximately triangular) profile is 50 degrees, or 100 bins. Thus, we estimate the total 
number of F/G Sgr stream stars in this 2.5 degree-wide stripe is between 30,000 and 50,000. 

In order to calculate the stream density at this position, we need to know the volume 
represented by the piece of the stream we have sampled. The FWHM of the stream is about 
25° in the observed stripe of sky. At a distance of 29 kpc from the Sun, this would imply 
a stream width of 13 kpc. [Note, however, that the Celestial Equator does not cut through 
the Sagittarius stream perpendicular to its path across the sky. The angle between the stream 
direction and the Celestial Equator is 28° |20j , thus implying that the actual width of the stream 
is (13 kpc) sin(28°) = 6 kpc] 

We estimate that the angle between the actual tidal stream and the projection of the tidal 
stream perpendicular to our line of sight is about 10° at this position in the sky so the 

distance along the stream in the plane of the Celestial Equator is (13 kpc)/ cos(10°) = 13 kpc. 
In the 2.5° wide stripe we sample 1.3 kpc perpendicular to the Celestial Equator. If the stream 
is in the range [4,8] kpc deep, then the volume is: 

V = (13 kpc)(1.3 kpc) ([4, 8]kpc) = [68, 140]kpc 3 . (6) 

Therefore, the density of F/G turnoff stars in that portion of the stream is: 

[30, 000; 50, 000] stars , „ 3 

Pf/G,82 = [68 ; i4Q]kpc J 3 = PIG, 740]stars/k P c 3 . (7) 



2.3.2 Estimating / 9 B hb,o//Obhb,82 

One now asks how much the stellar density varies as a function of position along the stream. 
This is not well known, but we will estimate it from BHB star counts pQ. It would be better 
if we could measure the number of F/G turnoff stars, as they a more stable indicator of stellar 
density. However, the only piece of the stream close enough to the Sun that this quantity has 
been measured reliably is the stripe 82 data discussed above. Table I summarizes the data and 
calculations described below. 

The stripe number (column 1) refers to the piece of sky scanned in the SDSS. Column two 
gives an estimate of the number of BHBs in the Sagittarius stream in each of those pieces of 



7 



data . The third column gives the distance from the Sun to the center of the stream, assuming 
the absolute magnitude for the BHB stars is M g = 0.7. 

We define 9 to be the angle between the normal to our line of sight towards a particular 
piece of the stream and the tangent to the stream at that point. We estimate 9 = 30°, 10°, 10° 
for stripes 10, 29, 82, respectively If = 0°, then the stream is perpendicular to our 

line of sight. If 9 is larger than that, then we see a larger distance along the stream in each 
solid angle of sky. The volume within the stream is inversely proportional to cos(#). cos(#) is 
tabulated in column 4. 

If we define <f> as the angle between the observations and the Sagittarius dwarf stream, 
then cos(</>) is found by dotting the Sagittarius dwarf plane normal vector with the SDSS 
observational plane normal vector. The normal to the plane of the Sagittarius dwarf tidal 
stream is: h = -0.064X + 0.970F + 0.233Z 0. In contrast to the referenced paper, we have 
defined X to be positive towards the Galactic center from the Sun. The stripes 10 and 82 are 
both on the Celestial Equator, so the normal to that observational plane is the north Celestial 
Pole: hpoie = -0.4834A" + 0.7472Y+ 0.4560Z. The normal to the plane of observation of stripe 
29 is: n 29 = -0.13091 + 0.4534F - 0.8816Z. If <j> = 90°, then we scan directly across the tidal 
stream. Smaller angles will net a larger effective volume and more stars; the volume is inversely 
proportional to sin(</>). We tabulate sin(0) = \/l — cos 2 (0) in column 5 of Table I. 

The column 6 in the table above gives a relative density of A stars in each of the three 
positions along the stream, corrected for projection effects (9), and the angle between the 
stream and the direction we scan across it (<j>). Following the calculation of F/G star density 
above, the volume sampled in each stripe, for a stream of cross sectional area A, is given by: 

v= ADtaM2Jn 
cos(0) sin 4> 

To estimate A star densities for a stream of diameter 6 kpc, divide the numbers in column 6 by 
tt(3 kpc) 2 tan(2.5°) = 1.23 kpc 3 . 

Comparing the numbers in column 6, we see that the density in the trailing tail is similar 
in the two places sampled. The density in the leading tail in the one place sampled is twice as 
high. The Sagittarius tidal debris at the position of the Sun is further along the leading tail, 
and we guess from these measurements that it could have a mass density similar to, or twice as 
large as, that measured in "stripe 82." 

Pbhb,q/pbhb,82 = [1,2]. (9) 
2.3.3 Estimating Af Sgr /A^ F/G , Sgr , M/L 

All that is left now is to find the conversion from the density of turnoff stars to a mass density 
of dark matter. We will do this by comparing the number of F/G stars in the Sgr dwarf 
galaxy to the estimated mass of the Sgr dwarf. Figure 6 of Newberg et al. {JO] shows a color- 
magnitude diagram of stars from two small fields near the center of the Sgr dwarf [22] , shifted 
into SDSS colors. The F/G stars in the stripe 82 data were selected with 0.0 < g — r < 0.6 
and 20.5 < g < 22.5. This color-magnitude selection was chosen to maximize the overdensity 
of Sagittarius stars in stripe 82. Since the F/G turnoff stars in Figure 6 of Newberg et al. [2*0] 
are a magnitude fainter in g, we selected all stars in that figure with 0.0 < g — r < 0.6 and 
21.5 < g < 23.5. There are 5786 F/G turnoff stars within these limits. 
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The field centers of the two Sgr dwarf fields 22 from which the data were taken were both 
within 1.5 degrees of the center of the Sagittarius dwarf. Figure 5 of shows the density of M 
giant stars along the major axis of the Sagittarius dwarf. Within 1.5 degrees of the center of the 
galaxy the density is reasonably flat. The surface brightness in V at the center of the Sagittarius 
dwarf is 25.2±0.3 magnitudes per square arcsecond |2IDE2]- The integrated apparent magnitude 
of the undisrupted portion of the Sgr dwarf is V = 3.63 Using a distance modulus of 16.9 to 
the dwarf, the total luminosity is 1.7 x 10 7 L©, comparable to previous estimates This 
means in one square arcsecond near the center of Sgr, the brightness is 4.25 x 10 8 times less 
than the brightness of the entire dwarf galaxy. So, if we found 5786 F/G turnoff stars in two 
holds 8.5' x 9' in size, then the total number of F/G turnoff stars in the Sagittarius dwarf is: 

The mass-to-light ratio in the Sgr tidal stream has not been measured by any technique. In 
analogy to other galaxies, we expect the luminous matter in the progenitor Sgr galaxy to be 
more concentrated than the dark matter. In tidal disruption, the outer portions of the dwarf 
galaxy are tidally stripped before the inner portions; the tidal radius decreases as the mass of the 
galaxy decreases. Therefore, one imagines that M/L in tidal streams is higher than M/L in the 
originating dwarf, and that the mass-to-light ratio increases with distance from the progenitor, 
as the more distant pieces were the first to be stripped. 

Figure 2 of Padmanabhan ct al. shows the radial profiles of luminous matter and dark 
matter for a stellar Hernquist profile embedded in an adiabatically compressed dark matter 
halo. M/L at the half light radius is twice the integrated M/L interior to that radius. At 1.5 
times the half light radius, M/L rises to three times M/L interior to the half light radius. At 
three times the half-light radius we virtually run out of luminous matter, so M/L becomes very 
large. 

Since our estimates of the luminous matter in the stream are of the same order of magnitude 
as the Sgr dwarf itself, the debris at the solar position probably was not stripped very far from 
the half light radius. We adopt as a reasonable upper limit that M/L in the stream could be a 
factor of three higher than the current M/L of the dwarf. 

Putting this all back into Eq. 0] we derive an estimated range for the local density of dark 
matter of: 

Ptail , = [210,740] x [1,2] x ^^of x [1,3] 
= [2, 190] x 10 4 A/ Q /kpc 3 

= [0.001, 0.07] GeV/cm 3 . (12) 

Comparing with the local density of dark matter not in the Sgr stream, — 0.3 GeV/cm 3 , 
we see that the Sgr stream contributes an additional 0.3%-23% to the total local dark matter 
density. 

Previously Stiff, Widrow, and Frieman |57] studied the tidal disruption of satellite galaxies 
falling into the halo of our own. They found that, with probability 0(1), the Sun should be 
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situated within a stream of density ~ 4% of the local Galactic halo density. Such a prediction 
is close to the density estimated for the Sgr stream if intermediate values for all variables are 
used. 



3 Direct Detection of WIMPs in the Stream 
3.1 Count Rates 

More than twenty collaborations worldwide are presently developing detectors designed to search 
for WIMPs. Although the experiments employ a variety of different methods, the basic idea 
underlying WIMP direct detection is straightforward: the experiments seek to measure the 
energy deposited when a WIMP interacts with a nucleus in the detector |28) . 

If a WIMP of mass m scatters elastically from a nucleus of mass M, it will deposit a recoil 
energy E — (p 2 v 2 /M)(l — cos 6*), where p = mM/(m + M) is the reduced mass, v is the speed of 
the WIMP relative to the nucleus, and 9 is the scattering angle in the center of mass frame. In 
this paper we consider the contribution to the detection rate from two dark matter components: 
the usual WIMPs in the Halo of our Galaxy plus the additional contribution from WIMPs in 
the Sagittarius tidal stream. We compute (following, e.g., 121311101) the differential detection 
rate per unit detector mass (i.e., counts/day/kg detector/keV recoil energy) as 

m = ^^[pWh(E,t)+p strV ME,t)] (13) 

where pt = 0.3 GeV/cm 3 is the standard local halo WIMP density (excluding the stream), 
Pstr — Ptaii,© is the dark matter density in the Sgr tidal stream, rjh is the mean inverse speed 
of WIMPs in the standard Galactic halo (excluding those in the stream), ?7 str is the mean 
inverse speed of WIMPs in the stream, er is the total nucleus- WIMP interaction cross section, 
q = v 2ME is the nucleus recoil momentum, and F{q) is a nuclear form factor that takes into 
account the loss of coherence in WIMP-nucleus interactions for momentum transfers comparable 
to or larger than the inverse nuclear radius (we normalize F(0) = 1). For detectors consisting 
of different nuclei, like e.g. Nal or CaW04, the total rate is a weighted average of the rates for 
the individual nuclei with weights equal to the mass fraction of each nucleus in the detector. 
We discuss these quantities in further detail in the following. 

The nuclear form factor F{q) depends on the type of WIMP-nucleus interactions, namely if 
they are spin-dependent or spin- independent, and reflects the mass and spin distributions inside 
the nucleus. For spin-independent interactions, the target nucleus can be approximated as a 
sphere of uniform density smoothed by a gaussian |31j , and the form factor follows as 

_ 3jsin(ggi) - qRi cos(qfli)] _ 9 v/2 n4 x 

*W~ q3R 3 ^> 

where Ri = (R 2 — 5s 2 ) 1 / 2 , s ~ 1 fm, and 

R ~ [0.91 (M/GcV) 1/3 + 0.3] x 10" 13 cm (15) 

is the radius of the nucleus. For spin-dependent interactions, the form factor is somewhat 
different but again F(Q) — 1. In general the form factor needs to be evaluated for specific 
detector nuclei. 
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For purely scalar interactions, 



4/< 



2 



Scalar = -!=- [Z f p + (A - Z)f n ] Z . (16) 



7T 



Here Z is the number of protons, A— Z is the number of neutrons, and f p and /„ are the WIMP 
couplings to nucleons. For purely spin-dependent interactions, 

a , sp i„ = (32A)Gf,/i 2 A 2 J(J + 1) . (17) 

Here J is the total angular momentum of the nucleus and A is determined by the expectation 
value of the spin content of the nucleus (see |28l El \TH I34j ) . 

For the estimates necessary in this paper, we take the WIMP-nucleon cross section to be 
cr p = 7.2 x 1CP 42 cm 2 , and take the total WIMP-nucleus cross section to be 4 

M .2 



(18) 

where the fi p is the proton- WIMP reduced mass, and A is the atomic mass of the target nucleus. 

Information about the WIMP velocity distribution is encoded into the mean inverse speed 
r,(E,t), 

rj(E,t)=[ ^A, (19) 

where 



ME , . 

*W = VV (20) 

represents the minimum WIMP velocity that can result in a recoil energy E and /a(u, t) is the 
(usually time-dependent) distribution of WIMP velocities u relative to the detector. 

In this paper, we neglect the rotation of the Earth and take the reference frame of the 
detector to be the same as that of the Earth. Thus we relate the velocity distribution /d(u, t) 
relative to the detector to the velocity distribution /(v) relative to the Galactic rest frame using 

/d(u,t) = /(v(t)), (21) 

v(i) = u + vg,(t), (22) 

where v®(i) is the velocity of the Earth in the Galactic rest frame. We discuss the motion of 
the Earth in the next section. 



3.2 Motion of the Earth in Galactic coordinates 

In this section we write an expression for the velocity of the Earth in the Galactic rest frame. We 
use a coordinate system in which X points toward the Galactic center, Y toward the direction 
of Galactic rotation, and Z toward the North Galactic Pole. 

For simplicity, we follow PU] and neglect the ellipticity of the Earth orbit and the non- 
uniform motion of the Sun in right ascension (an error of less than 2 days in the position of 

4 In most instances, f„ ~ f p so that the following equation results from Eq. 1161 . 



11 



the modulation maximum and minimum; more precise expressions for the motion of the Earth 
can be found in 35, 36 ). Thus we write the velocity of the Earth in terms of the Sun ecliptic 
longitude A(i) as 

v®(f) = v LSR + v Q + V® [ei sin X(t) - e 2 cos X(t)] . (23) 

Here 

v lsr = (0, 220, 0) km/s, (24) 

is the velocity of the Local Standard of Rest, which points in the direction of Galactic rotation, 
and 

v Q = (10, 13, 7) km/s (25) 

is the Sun's peculiar velocity. (The uncertainty in the Sun's peculiar velocity is of the order 
of 0.2 km/s in the Z direction and of as much as 3 km/s in the X and Y directions |37]; the 
corresponding uncertainty in the phase constant of the modulation is of the order of several 
days, larger than the error from neglecting the ellipticity of the Earth orbit.) Finally, 

V® = 27rA.U./yr = 29.8 km/s. (26) 

is the orbital speed of the Earth. 

The unit vectors ei and <32 define the plane of the Earth orbit, the ecliptic. They are in the 
direction of the Sun at the spring equinox and at the summer solstice, respectively. In Galactic 
coordinates, 

ei = (-0.0670, 0.4927, -0.8676), (27) 
e 2 = (-0.9931, -0.1170, 0.01032). (28) 

The Sun ecliptic longitude X(t) can be expressed as a function of time t in years with t — at 
January 1 as 

A(t) = 360° (t - 0.218). (29) 
Here 0.218 is the fraction of year before the spring equinox (March 21). 



3.3 Recoil Spectrum of Galactic and Sgr stream components 

In this paper we consider two different sources of WIMPs that contribute to count rates in 
detectors: WIMPs in the Milky Way halo and WIMPs in the Sgr stream. The function r](E, t) 
of Eas. (|13ll9|) is different for these two contributions. 

Galactic Component: For WIMPs in the Milky Way halo, the most frequently employed 
background velocity distribution is that of a simple isothermal sphere In such a model, the 
Galactic WIMP speeds with respect to the halo obey a Maxwellian distribution with a velocity 
dispersion ah truncated at the escape velocity v QSC , 

/,(v) = {ivd^ e ^' for l v ' < Vcsc (30) 
1 0, otherwise. 

Here 

v = y^2/3a h (31) 
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and 



N csc = erf(z) - 2zexp(-z 2 )/ 7 r 1 / 2 , 



with z — v esc /vo, is a normalization factor. For the sake of illustration, we take ah 
and v csc — 650 km/s. 

It can be shown that for the WIMPs in the isothermal halo 



(32) 
270 km/s 



Vh(E,t) 



1 



2v B (t)N a 
1 

2%(t)JV a 
0, 



erf (a; + y) 
erf (z) — erf (a; — y) 



erf (a; - j/) =ye~ 

\/7T 



for x < z 



:(z + y — x)e 



c , ( 33 ) 

tor z — y < x < z + y, 



for a; > 



where a; = u m in/wo, 2/ = v®(t)/vo, and z = v csc /v . 

Sgr Stream Component: For WIMPs in the Sgr stream, we can at first assume that they all 
move at the same velocity 



(0,0.233,-0.970) x 300km/s, 



(34) 



and later consider that in addition they have a small velocity dispersion er str somewhat larger 
than 20 km/s. The additional flux of WIMPs from the stream shows up as a 0.3-23% increase 
in the rate of nuclear recoils below the highest energy that WIMPs in the stream can impart 
to a target nucleus. Hence, there is a step in the energy recoil spectrum: the count rate in the 
detector is enhanced at low energies, but then returns to the normal value (due to Galactic Halo 
WIMPs) at all energies above a critical energy E c . This feature can be observed as a sharp 
decrease in the count rate at a characteristic energy that depends on the mass of the target 
nucleus, the mass of the WIMP, and the speed of the stream relative to the detector. 

If we neglect the velocity dispersion, the velocity distribution of the WIMPs in the Sgr 
stream becomes a Dirac delta function 



/(° r >(v) = £(v-v str ), 



(35) 



and their mean inverse speed computed from Eq. I|19(l is constant up to the characteristic energy 
E c (t), 



Here 9 is the Heaviside function, and 

U st r{t) = I 



0(E c (t) - E) 



Wstr(i) 



(36) 



(37) 



is the relative speed of the WIMPs with respect to the nucleus. The characteristic energy at 
which there is a step in the recoil spectrum is 



E c {t) 



2 ^ 2 r / + m2 



(38) 



This characteristic energy is the maximum recoil energy that can be imparted to the nucleus, 
and can be obtained as follows: The maximum momentum transferred from a WIMP to a 
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nucleus occurs when the WIMP bounces back and is q ma , x = 2/xii str (t). The maximum recoil 
energy of the nucleus then follows as E c {t) — q^ aax /(2M). 

Effect of Velocity Dispersion of Sgr stream: The effect of a velocity dispersion <7 str in the Sgr 
stream is to smooth out the edges of the step. We assume that the WIMPs in the Sgr stream 
follow a Maxwellian velocity distribution with bulk velocity v str , 



/str(v) 



1 



(27nr s 2 tl ./3) 3 /2 



-3|v-v Btr | 2 /2o 



(39) 



The mean inverse speed in this case can be obtained using the expression of rjh(E, t) after taking 
the limit t> csc — > oo and replacing v®(t) with u s tr(t) and vo with \j2ji cr str . 



r} stT {E,t) 



1 



2u st r(i) 



erf 



U s tr(t) 



erf 



u str (t) 



(40) 



The velocity dispersion of stars in the stream is observed to be roughly 20-30 km/sec [381 
1391 HU| . We will here take the velocity dispersion of the dark matter to lie in the same range, 
but caution that this number is very poorly known. 

3.4 Sgr Stream Annual Modulation 

It is well-known in the case of WIMPs in the Galactic Halo that the signal experiences an annual 
modulation ^E], with a peak on June 2 (for an isothermal halo). Similarly, the WIMPs in the 
Sgr stream also have an annual modulation in the count rate. As we will show, the count rate 
of the stream WIMPs peaks in June and is a minimum in December. In addition, the position 
of the step E c (t) depends on time, and so is annually modulated. It is higher in winter and 
lower in summer (i.e. its phase is opposite to that of the count rate). 
The time dependence of E c can be made explicit as follows. 



E c (t) = E^ {1 + A c cosM< - t c )]} . 



with u) = 2ir /l yr, 



E f) 
A c 

and t c is the solution of 

cos tot c ■ 

Here 



v 

M 



| Vstr - VLSR - V | 2 + Vi 



2VeV / af + a 



|v str - VLSR - V© | + V£ 



(1-2 



:, sincjt c = — 



a, 



e l • (Vgtr - vlsr - v ). 



(41) 

(42) 
(43) 

(44) 
(45) 



For our reference case of v stI — 300 km/s in the direction 0.233X — 0.970Z, we find that 
the amplitude of the modulation A c = 9.2% and that E c (t) is maximum on December 27 
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(t c = 0.991 yr) and minimum about six months later on June 28 [t c — 0.491 yr). Other values 
of w s tr between 100 km/s and 300 km/s, and other angles within 30° of the direction downward 
of the Galactic plane, give a modulation amplitude A c that can be as low as 2% or as high 
as 18%. The amplitude and phase of the E c modulation do not depend on the target nucleus, 
the WIMP mass, or the density of WIMPs in the stream. They depend on the magnitude and 
direction of the Sgr stream velocity. The direction of the stream will be known better and better 
as further observations of stars in the stream are made. 

The count rate of the signal from the stream in the detector is also modulated. For most of 
the recoil energy range, the phase of the count rate from the stream is opposite to the phase 
of the characteristic energy i.e. for our standard case it peaks on June 28. For energies near 
E c , there is a small range of energy recoils in which the count signal has the same phase as the 
characteristic energy and so peaks in December. The phase of the Sgr stream annual modulation 
is different from that of the normal annual modulation of WIMPs in the Galactic halo, which 
peaks on June 2. Contrary to the modulation of the characteristic energy, the count rate of 
the stream signal depends on the target nucleus, the WIMP mass, and the stream density in 
addition to the magnitude and direction of the Sgr stream velocity. 



3.5 Directional detection of WIMPs 

A possibility of separating events due to WIMPs in the Sgr stream from those in the local 
isothermal halo is to exploit detectors that can record the direction of individual nuclear recoils 
(directional detection). It has been difficult to build WIMP detectors sensitive to the direction 
of the nuclear recoils. However, a recent promising development is the DRIFT detector |41j . 
The DRIFT detector consists of a negative ion time projection chamber, the gas in the chamber 
serving both as WIMP target and as ionization medium for observing the nuclear recoil tracks. 
The direction of the nuclear recoil is obtained from the geometry and timing of the image of 
the recoil track on the chamber end-plates. Aim 3 prototype has been successfully tested, and 
a 10 m 3 detector is under consideration. The prototype contains ~ 0.2 kg of carbon disulfide 
(CS2) as active target. 

The differential detection rate including the direction of the nuclear recoil can be written 
using the results of 02| as 

- "' l [phA( w min,q) + Pstr^tr(«min,q)]) (46) 



dEdflq inm/i 2 

where q is a unit vector in the direction of the nuclear recoil, dd q is an infinitesimal solid angle 
in the direction q, v m i n is defined above fEa l20fl , and 

Tip, q) = / S(p - q • v) /(v) d 3 v, (47) 

is the 3-dimensional Radon transform of the velocity distribution function /(v) (in our case, of 
/fc(v) and /str(v), respectively). 

The Radon transforms of the velocity distributions we need have been computed in |42) . 
For the truncated Maxwellian distribution we assume for the WIMPs in the isothermal halo, 
we have 

MtW, q) = „ \ /2 _ [e-^+^W] 2 /^ _ e -vlJvl] . (48) 
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For the WIMPs in the Sgr stream, assuming a Maxwcllian velocity distribution with velocity 
dispersion a s t ri we have 

/strK^q) = ■L e- 3 ^^--'''^, (49) 
V 27T<J h/ s 

where 

U str = V str - V©(t). (50) 

When the velocity dispersion is neglected, the Radon transform for the stream reduces to 

/rtr (^min, q) = ${Vm.in - 5 ' U str ) (51) 

The latter case has a simple geometrical interpretation in the space of vectors v m [ n q. The 
Dirac delta function forces the velocity w m ; n to change with the recoil direction q according 
to fmin = q • u str . As q varies, the vectors u m i n q describe the surface of a sphere centered in 
2 u s tr and passing through the origin. The sphere has diameter |u str |, and it is symmetric for 
rotations around u str - 

With a small velocity dispersion in the Sgr stream, the recoil vectors iWnq due to WIMPs 
in the stream are still approximately concentrated around the sphere just described. On the 
contrary, the recoils due to WIMPs in the isothermal halo are smoothly distributed over a wide 
region of w m i n q space. It is this big difference in the recoil distributions from WIMPs in the 
isothermal halo and in the Sgr stream that may allow their separation. 

For detectors consisting of different nuclei, like e.g. CS 2 for DRIFT, the total rate is a 
weighted average of the rates for the individual nuclei with weights equal to the mass fraction 
of each nucleus in the detector. In this case it is not possible to know with which type of 
nucleus the WIMP scattered in a given event. Hence it is not possible to obtain the magnitude 
of the nuclear momentum q — \J2ME or the velocity w m i n = J ME/ (2/i 2 ), since they depend 
on the nuclear mass M. In this case, it is convenient to present the directional recoil spectrum 
dR/dEdVL q in terms of the quantity 

E = Eq, (52) 

which is a vector with magnitude equal to the recoil energy and direction pointing in the recoil 
direction. We call it the directional recoil energy. The (X, Y, Z) components of the directional 
recoil energy are the product of the recoil energy and the direction cosines of the nuclear recoil 
momentum. Since the volume element in the E space is E 2 dEdfl q , the relevant rate to plot is 

J—**- (53) 

e 2 dEdn q ' { ' 

which is measured in counts per kg-day per keV 3 . 



4 Results 

In Table II, we list the values of the characteristic energy E c of the step in the energy recoil 
spectrum due to WIMPs in the Sgr stream appropriate to our reference case (v s t r = 300 km/sec 
toward (l,b) = (90°,— 76°)) for different target nuclei and WIMP masses. As discussed above, 
E c experiences an annual modulation, with a maximum in June and a minimum in December. 
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Target nucleus 


m=60 GcV 


m=100 GeV 


m=500 GeV 


£; c (Dcc27) 


£ c (Jun28) 


£ c (Dec27) 


£ c (Jun28) 


£ c (Dec27) 


£ c (Jun28) 


V3 Ge 


42.3 


35.1 


68.2 


56.7 


149.7 


124.3 


2a Na 


33.0 


27.4 


41.2 


34.3 


56.0 


46.5 




37.8 


31.4 


70.2 


58.3 


219.6 


182.4 


4U Ca 


40.2 


33.4 


56.1 


46.6 


91.7 


76.2 


183 W 


32.6 


27.1 


65.9 


54.7 


267.0 


223.4 


16 


27.1 


22.5 


32.1 


26.6 


40.0 


33.2 


i3i Xe 


37.5 


31.1 


70.0 


58.2 


223.8 


185.9 




31.6 


29.4 


45.6 


37.9 


64.8 


53.8 



Table 2: Characteristic recoil energies E c in keV from WIMPs in the Sagittarius stream. Note 
that the value of E c has an annual modulation with a peak on December 27 and a minimum 
on June 28. Numbers have been obtained for our reference case, with a stream velocity of 300 
km/sec toward (l,b) = (90°, -76°). 



The annual modulation of the characteristic energy E c can in principle be used to identify 
the signal due to WIMPs in the Sgr stream. In particular, the annual modulation is a useful 
tool to differentiate between an increase in the count rate due to WIMPs in the Sgr stream at 
energies below E c , and an increase due to modifications of the velocity distribution of Galactic 
halo WIMPs. These possibilities are difficult to distinguish by simply examining the shape of 
the recoil spectrum because of the finite energy resolution and the binning of data in recoil 
energy in a real experiment. For this tool to be effective, the energy resolution must of course 
be better than the difference between the maximum and minimum values of E c . As discussed 
in Section 3.4 above, this difference corresponds to an energy resolution of 2-18% depending on 
the exact velocity of the Sgr stream. For our reference case of v s tr = 300 km/sec in the direction 
(I, b) = (90°, —76°), the value of E c is modulated by 9% over the course of the year. As long 
as the uncertainty in energy due to the energy resolution is smaller than the energy shift in E c , 
one should be able to detect its annual modulation. 

We note that the values of E c are indeed in the range accessible to current and upcoming 
dark matter experiments. For example, the Nal detector of the DAMA experiment 031 nas an 
energy threshold of 2 keV electron equivalent, which corresponds to 22 keV recoil energy for 
iodine; since iodine is much heavier than sodium, WIMP interactions with iodine dominate the 
count rate, and the conversion factor between keVee and keV is a factor of 11 for iodine. Hence 
for WIMP masses heavier than 50 GeV, the step lies above the DAMA threshold energy, and 
DAMA can in principle detect it. The question is whether or not DAMA can in fact identify 
an increase in count rate with a signal from the Sgr stream. 

In Figure 1, we have plotted the count rate of 60 GeV WIMPs in a Nal detector such as 
DAMA as a function of recoil energy. The dotted lines (towards the left) indicate the count rate 
due to Galactic halo WIMPs alone for an isothermal halo. The solid and dashed lines indicate 
the count rate including the WIMPs in the Sgr stream for our reference case of v str = 300 
km/sec in the direction (l,b) = (90°,— 76°) with a stream velocity dispersion of 20 km/sec. 
In this plot, we have taken the Sgr stream to contribute an additional 20% to the local Halo 
density of our Galaxy. As discussed earlier, the best estimate of the contribution varies from 
0.3% to 23%. The solid and dashed lines are for June 28 and December 27 respectively, the 
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2 4 6 

E ee [keV] 

Figure 1: Count rate of 60 GeV WIMPs in a Nal detector such as DAM A as a function of recoil 
energy. The dotted lines (towards the left) indicate the count rate due to Galactic halo WIMPs 
alone for an isothermal halo. The solid and dashed lines indicate the step in the count rate 
that arises if we include the WIMPs in the Sgr stream for v s t r = 300 km/sec in the direction 
(l,b) — (90°,— 76°) with a stream velocity dispersion of 20 km/sec. The plot assumes that 
the Sgr stream contributes an additional 20% of the local Galactic halo density. The solid 
and dashed lines are for June 28 and December 27 respectively, the dates at which the annual 
modulation of the stream is maximized and minimized. 
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dates at which the stream annual modulation is maximized and minimized. As noted above, 
the annual modulation of the value of E c for this reference case is a 9.2% effect, and peaks at 
a different date than the standard annual modulation in the Galactic halo (which causes the 
count rate to peak on June 2). 

The DAMA experiment has an enormous exposure of 107731 kg-days. For a 60 GeV WIMP 
mass and a maximal Sgr stream that is 20% of the local halo density, we compute the following 
number of events expected in 1 keV bins (2-3 keVee, 3-4 keVee, 4-5 keVee, 5-6 keVee), averaged 
over a year: (136744, 51814, 20873, and 9367) with the Sgr stream, and (98370, 45745, 20717, 
and 9214) without the Sgr stream. 5 The location of the step is roughly at 3 keVee. For a 60 
GeV WIMP mass and a stream that is 20% of the local halo, in the 2-3 keVee bin, the difference 
in number of counts with and without the stream is 38374, which is 104 times the square root 
of the number of counts with the stream. Hence we find that, in the 2-3 keVee bin, the presence 
of the stream is detectable at the 104<r level. In the 3-4 keVee bin, we find that the stream 
is detectable at the 27<r level. These results have assumed that the density of WIMPs in the 
Sgr stream is 20% of the local halo. If, as is more reasonable, the density of WIMPs in the 
stream is instead 4% of the local halo, the stream is detectable at the 24<r and 6<r level in the 
2-3 and 3-4 keVee bins respectively. In more detail, if the stream density is 4% of the local halo, 
the number of expected events in 1 keV bins from 2 to 6 keVee, averaged over a year with an 
exposure of 107731 kg-days, is (106045, 46959, 20748, and 9245). The expected count rates for 
other stream densities can easily be obtained by scaling the stream contribution in the numbers 
quoted above, for the contribution of the stream scales proportionally with the ratio of stream 
density to local halo density. 

We note that, in the presence of a stream, the date of the maximum count rate changes. For 
the reference case of v s t r = 300 km/sec in the direction (I, b) = (90°, —76°) the count rate due to 
the WIMPs in the stream alone peaks on June 28 and December 27 at different values of recoil 
energy. Depending on the contribution of stream WIMPs relative to the local halo WIMPs, 
the overall count rate may peak between these dates and the usual June 2. For a Sgr WIMP 
density that is 20% of the local halo, the overall count rate peaks on June 28 and December 27 
for different values of the recoil energy. For lower values of the stream density, the date of peak 
count rate changes and eventually returns to the same value as without a stream. A further 
investigation of this sensitivity is in progress. 

The sapphire (AI2O3) and calcium tungstate (CaW04) detectors of CRESST have a 
better energy resolution. Its smaller detector mass requires a longer exposure. The energy 
threshold is low, ~ 1 keV. 

The 73 Ge detector of the CDMS experiment has a low energy threshold at 10 keV and a 
good energy resolution, about 0.5 keV at threshold and about 1.5 keV at 100 keV. As can be seen 
from Table II and Figure 1, which shows the energy recoil spectrum in 73 Ge, the characteristic 
energy of the reference model varies from 25-120 keV, easily above the energy threshold of the 
experiment. In addition, the CDMS energy resolution in this range varies between 2% and 
4%, so the CDMS experiment has the required energy resolution to be sensitive to the annual 
modulation of the Sgr stream for most of the range of expected velocities. The EDELWEISS 
experiment also with a 73 Ge detector, has an energy threshold of 20 keV. 

5 For these stream parameters, iodine does not provide a signal from the stream beyond 5 keVee, while sodium 
gives rise to a long but weak tail of stream events extending up to 21 keVee. This tail, however, contributes 
less than one event per keVee in 107731 kg-days beyond 11 keVee. The long tail is due to the very different 
quenching factors of iodine and sodium: 0.09 and 0.3, respectively. 
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CDMS-II at Soudan is expected to have an exposure of 2500 kg-day 07] ■ Using the same 
parameters as in Figure 2, in particular a WIMP mass m — 60 GeV, we find the following 
number of events in 10 keV bins (10-20 keV, 20-30 keV, 30-40 keV, 40-50 keV, 50-60 keV), 
averaged over a year: (3083, 2138, 1386, 666, 383) with the Sgr stream and (2473, 1592, 1006, 
623, 383) without the Sgr stream. The last bin has the same number of counts with or without 
the stream because it is well above the characteristic energy of the step (in this case, E c varies 
between 35.1 and 42.3 keV). In the first bin, the difference in number of counts with and without 
the stream is 610, which is 11 times the square root of the number of counts with the stream. 
Hence we find that, in the first bin, the presence of the stream is detectable at the 11<t level. 
In the second and third bins we find that the stream is detectable at the 12cr and lOcr level, 
respectively. This assumed that the density of WIMPs in the Sgr stream is 20% of the local 
halo density. With 4% density, instead, the stream should be detectable in the same bins at 
the 2.4a, 2.7a, and 2.3<r levels, respectively. With the proposed CryoArray detector, which will 
have 500000 kg-day of data @7J, the Sgr stream should be easy to identify. 

In Figure 2, we have plotted the count rate of 60 GeV WIMPs in a 73 Ge detector (such 
as CDMS or Edelweiss) as a function of recoil energy. The parameters of the curves are the 
same as in Figure 1. The dotted lines (towards the left) indicate the count rate due to Galactic 
halo WIMPs alone for an isothermal halo. The solid and dashed lines indicate the count rate 
including the WIMPs in the Sgr stream for our reference case. The characteristic energy of the 
step is 35.1 keV in June and 42.3 keV in December. 

Upcoming detectors XENON ggj and ZEPLIN-II gg] have low energy thresholds of 2-4 keV. 
With 1 ton detectors made of liquid xenon planned, these experiments should see enormous 
count rates and be able to distinguish the Sgr stream. 

In Figure 3 we illustrate the power of directional detectors. The figure shows the count rate 
of 60 GeV WIMPs in a CS2 detector (DRIFT) as a function of recoil energy and direction of the 
nuclear recoil. As above, we have assumed that the Sgr stream impinges on the Galactic plane 
with a speed of 300 km/s in the direction (I, b) = (90°, —76°). We plot the rate at the maximum 
of the stream modulation on December 27. The figure shows the count rate in a 2-dimensional 
slice of the 3-dimensional recoil space. The chosen slice is perpendicular to the direction of 
Galactic rotation and defined by a recoil energy of 20 keV in that direction. The horizontal axis 
represents recoils in the direction of the Galactic center (left) and Galactic anticenter (right); the 
vertical axis represents recoils in the direction of the North Galactic Pole (upward) and South 
Galactic Pole (downward) . The gray scale indicates the count rate per kilogram of detector per 
day and per unit cell in the 3-dimensional energy space. Lighter regions correspond to higher 
count rates. The white band on the upper part is the location of nuclear recoil due to WIMPs 
in the Sgr stream. The fuzzy gray cloud at the center contains recoils due to WIMPs in the 
local isothermal halo. The two WIMP populations can in principle be easily separated, given 
a sufficient exposure. Figure 4 shows a section of Figure 3 (as defined in the captions) and 
illustrates that the Sgr stream clearly stands out above the Galactic halo WIMPs. 

We can estimate the exposure needed to identify the Sgr stream by considering binning 
directional data into cells in directional energy space of size 10 keV. The count rate per bin is 
approximately given by the rate in figures 3 and 4 multiplied by the volume of the cell (10 3 
keV 3 ). The count rate for the Sgr stream results in ~ 10~ 2 counts/kg-day/cell, and the count 
rate for the isothermal component is ~ 3 x 10~ 3 counts/kg-day/cell. Using a \HV argument, we 
see that the Sgr stream may be distinguished at the 3cr level with an exposure of ~ 10 4 kg-days. 
For DRIFT, this corresponds to exposing a 30 m 3 of CS2 for 5 yr. This is in the upper range 
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Figure 2: Count rate of 60 GeV WIMPs in a 73 Ge detector as a function of recoil energy. The 
dotted lines (towards the left) indicate the count rate due to Galactic halo WIMPs alone for an 
isothermal halo. The solid and dashed lines indicate the step in the count rate that arises if we 
include the WIMPs in the Sgr stream. The parameters of the curves are the same as described 
in Figure 1. 
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Figure 3: Count rate of 60 GeV WIMPs in a CS 2 detector (DRIFT) as a function of recoil 
energy and direction of the nuclear recoil. The figure shows the count rate in a 2-dimensional 
slice of the 3-dimcnsional recoil space. The chosen slice is perpendicular to the direction of 
Galactic rotation and defined by a recoil energy of 20 keV in that direction. The horizontal axis 
represents recoils in the direction of the Galactic center (left) and Galactic anticenter (right); the 
vertical axis represents recoils in the direction of the North Galactic Pole (upward) and South 
Galactic Pole (downward) . The gray scale indicates the count rate per kilogram of detector per 
day and per unit cell in the 3-dimensional energy space. Lighter regions correspond to higher 
count rates. The white band on the upper part is the location of nuclear recoil due to WIMPs 
in the Sgr stream. The fuzzy gray cloud at the center contains recoils due to WIMPs in the 
local isothermal Galactic halo. The two WIMP populations can in principle be easily separated, 
given a sufficient exposure. 
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Figure 4: Section of figure 3 with E(\x = 0. 
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of the current plans for DRIFT-2, whose volume is expected to be in the range 10-30 m 3 . 
Therefore DRIFT has the capability of identifying WIMPs in the Sgr stream. 



5 Conclusions 

Recent observations of the Sagittarius dwarf spheroidal galaxy (Sgr) indicate the existence of 
tidal streams of stars that pass through the solar neighborhood. If the mass-to-light ratio in the 
streams is at least comparable to that of the main body of the Sgr galaxy, then one can detect 
the existence of dark matter in the streams. Under the assumption that the stream passes 
through the solar position, the dark matter stream plausibly has a density ~(0.3-23)% of the 
local density of our Galactic Halo. If dark matter consists of WIMPs, the extra contribution 
from the stream gives rise to a step-like feature in the energy recoil spectrum in direct dark 
matter detection, i.e., an increase in the rate in detectors at energy recoils less than some critical 
value E c that depends on the target nucleus and on the time of year. The location of the step 
experiences an annual modulation that will be useful in identifying the existence of the stream. 
The count rate in the detector is also modulated annually, with a maximum on June 28 and 
minimum on December 27 except near the characteristic energy where the phase is opposite. 

In the CDMS experiment, for our best estimate of stream velocity (300 km/sec) and direction 
[(/,&) — (90°,— 76°)], the location of the step oscillates yearly between 25 and 35 keV, with a 
phase opposite to that of the count rate. With two kg of Germanium, the stream should be 
detectable at three er's in four years of data with 10 keV energy bins. Planned large detectors 
like XENON, CryoArray and the directional detector DRIFT may be able to identify the Sgr 
stream. 
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